Electron induced light emission in photonic crystals 

L. A. Blanco^'* and F. J. Garcia de Abajo^'^ 
^Donostia International Physics Center (DIPC), 
Aptdo. 1072, 20080 San Sebastian, Spain 

^Centra Mixto CSIC-UPV/EHU, 
Aptdo. 1072, 20080 San Sebastian, Spain 
(Dated: February 1, 2008) 

Abstract 
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The interaction of a fast electron with a photonic crystal is studied by solving the Maxwell equations 
exactly for the external field provided by the electron in the presence of the crystal. The polarization 
currents and charges produced by the passage of the electron give rise to the emission of the so-called 
Smith-Purcell radiation. The emitted light probability is obtained by integrating the Foynting 
vector over planes parallel to the crystal at a large distance from the latter. Both reflected and 
transmitted light components are analyzed and related to the photonic band structure of the crystal. 
Emission spectra are compared with the energy loss probability and also with the reflectance and 
transmittance of the crystal. 
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I. INTRODUCTION 



Photonic crystals have received considerable attention as many technological applications 
are being proposed including control of atomic and molecular spontaneous emission or 
light guiding and confinement [0, ^, |, |[. The main physical aspect that is used for these 
purposes is the existence of photonic bandgaps in such systems, that can lie in the near 
infrared ||^, ^ ^ or in the visible , which allow light to follow a determined path depending 
on the structure of the crystal. Despite recent advances in the production of photonic crystals 
that operate in both regions of the light spectrum , fabrication defects can still limit their 
applicability. However, a quantitative determination of the number and type of the defects 
in a crystal is difficult to perform. Electrons seem to be a valid probe to that end. Here, we 
study the interaction of a fast electron with a crystal as a first approach to this problem. 

In particular, a fast electron moving near a photonic crystal can polarize its constituents 
inducing charges and currents that cannot follow the electron motion and suffer acceleration 
as they evolve in the inhomogeneous crystal, giving rise to light emission. This is similar to 
the so-called Smith-Purcell effect fll], |12|, [13], 0, [T|, [l^] . 

The present work is intended to provide a new way to characterize photonic crystals 
by observing the spectrum of the light emitted when electrons are travelling parallel to the 
surface of the crystal. The electric field of the moving electron is decomposed into evanescent 
plane waves |T8|, and the outgoing light is produced by the diffraction of these waves in the 
crystal. This procedure is sketched in Section II. Numerical results for a complete-band-gap 
photonic crystal formed by air spheres in silicon (inverted opal) are shown in Section HI. 
We use electrons travelling within air and silicon, in order to reflect different electrodynamic 
effects, although a possible experimental observation could be more difficult in the latter 
case. Gaussian atomic units will be used from now on, unless otherwise specified. 



II. THEORETICAL FRAMEWORK 



A theoretical description of the interaction between an electron and a crystal is given next. 
The crystal will be composed of several layers perpendicular to the z direction, beginning 
in 2; = and extending towards the negative-z region. Outside the crystal a medium with a 
dielectric function e will also be considered. All media will be assumed to be non-magnetic 
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(/i = l). 

The electron will be considered to follow a trajectory described by rj = {vt,yo,ZQ), with 
Zq > 0, and its electric field is found to be, in frequency space lj and in the absence of the 
crystal, 

Eo(r, u) = J dQy eiK±-[r-(o,yo,.o)])E± , (1) 

where Kq = (Q,±irQ), Q = {uj/v,Qy) is the two-dimensional momentum parallel to the 
surface, and Tq = Q^ — uJ^ejc? , with Rejrq} > 0. In this expression, Eq exp(iKQ ■ r) is 
a plane wave, which can be expressed as a sum of s and "p components, Eq = s Gq ^ + 
i?Q p Bq p. The + (— ) sign in these expressions stands for a wave moving towards the positive- 
(negative-)z region. When the electron moves in vacuum, all of these waves are evanescent 
(Fq is real and positive). However, when it moves in a medium described by e > and 
V\fe > c, some of those waves represent propagating Cherenkov radiation. 

Due to the crystal symmetry, an incident wave with momentum Kq will only produce 
a discrete set of reflected (transmitted) waves of momentum Kq_^q (Kq^q), G being a 
reciprocal surface lattice vector. Therefore, the reflected and transmitted wave amplitudes 
can be expressed as 

a' 

[^Q+G,J* = H^QgI^Q,^']*) (2) 

where a and a' run over polarizations s and p, the super-indices r, t and i stand for reflected, 
transmitted and incident components, respectively, and -Rqg and Tqq are the reflection and 
transmission coefficients of the crystal. These coefficients are calculated using the layer KKR 
method |T^, in which the self-consistent electric field is constructed in terms of multipole 



expansions around the crystal objects by multiple scattering within each layer and later 
among layers. 

Each of the incident plane waves in Eq. ([l|) gives rise to a set of reflected and transmitted 
waves whose amplitudes are obtained according to Eq. (H). The electric fleld in the positive- 2; 
region can be constructed as the sum of Eq and the reflected fleld, 

E(r, uj) = Eo(r, u) + E^(r, u), 
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where 

G.o-cr' 

A similar expression is found for the transmitted electric field at the other side of the crystal 
(negative z's): 

G.crcr' 

The light emission probability per unit length is calculated by means of the projection 
of the Poynting vector V over the z and —z directions (for reflected and transmitted com- 
ponents, respectively), which is integrated over the time and over a plane parallel to the 
crystal surface. The resulting expression can be written 

J dtj dxj dyV ■ {±z) = ujdujP{uj), (3) 

where P{uj) can be interpreted as the emission probability per unit length and per energy 
range uj. The average over parallel impact parameters has also been performed. When the 
Cherenkov condition v'^e > is not fulfilled, the direct waves of Eq. (|l]) are all evanescent, 
and the light emission by reflection and transmission at the crystal (i.e. emission towards the 
2; > and 2; < regions, respectively denoted reflected emission and transmitted emission 
hereafter) reduces to 
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where 



_ iQyUJ k 

Aq = irql 

and the prime in the summation over G indicates that only non-evanescent outgoing waves 
(Tq+g < 0) Kiust be included. 



On the other hand, if f^e > c^, there exists an interference between the direct Cherenkov 
radiation and the reflected field. The emission probability expressions take a more involved 
form in this case: 



27rA;- 



a G 

dQy AQ|/.|2 + 2Re(AQ^rQo<,e-2^^'^ 
+ / rfQ,EVQG.l'^Q+Ge-'^'^^« , 

•IQo G J 



E' 

G 



^ QGcrpAq+G 



+ 



27ik- 



E' 

o-G 



y I ^ I 



2/FQGcr| ^Q+gC 



+ 



Z' |2 
+ / dQy\tQGa\ Aq+g + 

-'-Qo 

+ r dQy\tQGa\^Aci+Ge-^''^''] 

•'Qo J 



IQo 

where Qo = ^t — uP'lv'^, fp = —k/ (vQy/e) and fs = iQyUJ j {QVc^c^). The integrands in these 
expressions describe the reflected and transmitted wave components for external evanescent 
{\Qy\ > Qo) and propagating {\Qy\ < Qo) waves. In particular, the first term inside the 
square bracket of Pr reduces to the Cherenkov emission probability in a bulk homogeneous 
medium, while the remaining terms describe reflected and interference components. 



III. RESULTS AND DISCUSSION 



We have applied the previous formalism to the case of a 100-keV electron moving parallel 
to a slab consisting on several fcc(lOO) and fcc(lll) planes of air spheres in silicon. The gap 
characteristics of this system have been theoretically described in : a complete photonic 
bandgap is opened when the ratio between the sphere radius r and the cubic lattice constant 
a lies in between 0.335 and 0.374. We have chosen an intermediate value of r/a = 0.342, 
and the frequency range is taken close to the communications wavelength, A = 1.55/im, in 
which case esi = 11.9. 

Results for crystals composed of 8 fcc(lll) planes are shown in what follows, although 
we have also carried out calculations for 8 fee (100) planes and the results are qualitatively 
similar. Two different types of calculations are presented, a) In Fig. 1, the crystal is 
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embedded in Si, so that the electron moves within sihcon outside the crystaL In this case, 
the electron produces Cherenkov radiation even without the presence of the crystal whenever 
the electron velocity is larger than the speed of light in the medium, which is the case at 
100 keV. b) In Fig. 2, the crystal is embedded in air and the electron moves also outside 
the crystal, that is in air. Therefore, no direct Cherenkov radiation is produced. 

As the total energy must be conserved in the present case in which the material is trans- 
parent (that is, e is real) the sum of the reflected and transmitted probabilities must be the 
same as the electron energy loss probability, which has been calculated from the retarding 



force exerted on the electron by the induced electric field, as discussed elsewhere ^ 



The main difference between the two cases considered in Figs. 1 (a) and 2 is that when 
Cherenkov radiation is emitted the energy loss and the reflected emission are very much 
enhanced [Fig. 1 (a)]. This is due to the fact that the Cherenkov emission is dominant as 
compared to the diffraction of evanescent waves. For the transmitted light, as expected, the 
emission probability within the gap frequencies is strongly suppressed in both cases. 

When the electron travels in air (Fig. 2), the energy loss and the reflected light proba- 
bilities are much smaller than in the previous case. However, the order of magnitude of the 
transmitted field is the same as in the Cherenkov case. Notice that huge peaks are produced 
in the reflection curve near the bandgap region, probably connected with band edge effects. 
This type of energy losses are well below the current energy resolution of STEM machines 
(around 0.5 eV), but since the sample is transparent, any energy loss must be converted into 
photon emission, for which the resolution is much better, as shown in and sufficient to 
resolve the features discussed in this work. The integrated area of the sharp peaks in Fig. 2 
amounts to approximately 5 x 10~^ photons per electron per nm, which for path lengths of 
the order of a several microns, would result in a measurable photon intensity of a percentage 
of the number of electrons. 

In both cases the gap regions agree reasonably well with those found for the transmission 
of external electromagnetic waves under normal and Cherenkov-angle incidence conditions 
as well as with the band structure of the system, which are shown for the first case (electron 
moving in Si) in Figs. 1 (b)-(c). Note that the transmitted emission gap [Fig. 1 (a)] is wider 
than the complete gap (vertical dashed lines throughout all the Fig. 1), since the latter 
corresponds to a particular selection of incident directions (the Cherenkov cone). The same 
is true for the transmittance [Fig. 1(b)]; actually, the case of normal incidence corresponds 
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to the r point. The depletion of hght transmission within the gap could be used to detect 
defaults (stacking faults, missing atoms, etc.), whose relatively small contribution to the 
emission would be amplified because that would be the primary origin of light transmitted 
within the gap. 

Finally, in Fig. 3 we show the dependence of the light emission probability on the distance 
between the electron and the crystal, for a given frequency within the bandgap. As expected, 
both the transmitted and reflected probabilities decrease when increasing the electron-crystal 
separation, but it is interesting how the magnitude of the interference between the direct 
(Cherenkov) and the reflected flelds produces an oscillatory pattern at the front side of the 
crystal. For the transmitted field, the probability reaches an asymptotic value which comes 
from the transmission of only Cherenkov incident waves. 
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FIG. 1: (a) Frequency variation of the energy loss probability (thick solid curve), the 
reflected plus direct light emission probability (dashed curve) and the transmitted light 
emission probability (thin solid curve) for a 100-keV electron moving in front of 8 fcc(lll) 
planes of air spherical voids of radius r=418 nm in silicon (e = 11.9), where the lattice 
constant a is 1222 nm (the filling fraction is 67% and the separation between sphere centers 
is o?=873 nm). The electron is moving inside Si at a distance of 876 nm from the void spheres, 
(b) Transmittance for hght inciding normally (thick curve) and with the Cherenkov angle 
9c — 68.1° (thin curve) for the same system, (c) Bulk band structure projected on the (111) 
plane for the same system, where the shaded zones indicate regions where no propagating 
electromagnetic waves are allowed. A full gap is observed in the range uod/c ~ 3.23 — 3.35. 

FIG. 2: Frequency variation of the energy loss probability (thick solid curve), the reflected 
plus direct light emission probability (dashed curve) and the transmitted light emission 
probability (thin sohd curve) for a 100- keV electron moving in front of 8 fcc(lll) planes of 
air spherical voids of radius r=418 nm in silicon (e = 11-9), where the lattice constant a 
is 1222 nm (the fiUing fraction is 67% and the separation between sphere centers is d=873 
nm). The surface of the crystal is 131 nm from the last sphere. The electron is moving in 
air at a distance of 800 nm from this surface. 

FIG. 3: Variation of the reflected plus direct emission probability (dashed curve) and the 
transmitted emission probability (solid curve) with the distance between a 100-keV electron 
and the first of the 8 fcc(lll) planes of air spherical voids of radius r=418 nm in silicon, 
where the lattice constant a is 1222 nm. The emitted hght wavelength is A = 1.68//m (i.e, 
ud/c — 3.25). The asterisk (*) corresponds to the electron-spheres distance of Fig. 1 (a). 
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